Mean escape time over a fluctuating barrier 
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An approximate method for studying activation over a fluctuating barrier of potential is proposed. 
It involves considering separately the slow and fast components of barrier fluctuations, and it applies 
for any value of their correlation time r. It gives exact results for the limiting values r and 
r ^ oo, and the agreement with numerics in between is also excellent, both for dichotomic and 
Gaussian barrier perturbations. 
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Ever since Kramers seminal paper [3| the fluctuational 
escape over a potential barrier has been a paradigm for 
a thermal activation process. Recently, activation in the 
presence of time- varying fields have become a subject of 
great interest due to the discovery of many counterin- 
tuitive noise-assisted effects, like stochastic resonance '5| 
or transport in Brownian motors The nonequilib- 
rium character of these problems hinders, however, the 
direct application of many ideas and methods developed 
for investigation of the static Kramers problem 0| (e.g., 
detailed balance or rate concept). On the other hand, 
as the time-scale of variation of the driving signal is in- 
dependent of the internal dynamics of the system, stan- 
dard adiabatic methods are restricted to certain ranges of 
parameters, only. Hence, an approach which overcomes 
these difficulties and applies for the whole range of time 
variability of the perturbation, is of great importance. 

In this letter we address this problem for an activation 
over a randomly fluctuating barrier. The subject is inter- 
esting not only due to its ubiquity in many branches of 
physics, e.g., in relation to ligand binding to heme pro- 
teins IeJ], transport processes in glasses or dye laser 
with a fluctuating pump parameter Q , but especially be- 
cause of the phenomenon of resonant activation Q — the 
appearance of a minimum of the mean activation time T 
as a function of the correlation time t of barrier fluctu- 
ations. The dependence T{t) can be calculated exactly 
merely for simple models l a lSlijlCrl , for more general cases 
the approaches 0, 0, Il3lllini5.j proposed till now ap- 
ply to some ranges of r, only. Irrespective of technical 
differences they are all based on the rate concept, which 
assumes a quasi-stationary equilibrium before the activa- 
tion happens and applies for r ^ T, and/or kinetic de- 
scription for r ^ \n{l/q) {q states for the thermal noise 
intensity) when the escape events are uncorrelated with 
the potential variations. Although for small enough q 
in an extended region hi{l/q) <ti t <^ T both approx- 
imations coexist and give similar results neverthe- 
less the proper smooth connection between them remains 
the main theoretical challenge. Below, we present an ap- 
proach which is valid for any r. It gives exact values of 
T(t) in the limits t and r — s- oo, and a very good 
approximation in between. 



We study an overdamped Brownian particle driven 
by a (thermal) Gaussian white noise ^(t) of zero mean, 
which moves in a stochastically varying potential. Its 
static part U{x) has a monostable or bistable form and 
the random part V{x)z{t) is generated by a station- 
ary Markovian noise z{t) of zero mean and correlation 
C(t) = Q/rexp(-|t|/T). Following [ll[l3 we assume a 
general form for its intensity Q(r) = Qqt" {0 < Qo ~ 
const , < a < 1), which gives the mostly studied cases 
with r-independent intensity (a = 0) or variance (a = 1) 
as special cases. Two types of z{t) are considered: an 
Ornstein-Uhlenbeck noise (OUN), which is Gaussian with 
variance D — Q/t, and a dichotomic noise (DN), which 
flips between two values ±V1d with the rate 7 = 1/(2t). 
Although they essentially differ — the former is continu- 
ous, the latter discrete — nevertheless, they influence the 
activation process very similarly |lj and the main steps 
of the presented description are the same. The dynam- 
ics of the system is given by the nonmarkovian Langevin 
equation 



V'{x)zit)+m 



(1) 



Introducing the two-dimensional markovian stochastic 
process {x{t), z{t)} one can formulate the evolution equa- 
tion for the joint probability distribution P{x, z, t): 



di 



P{x, z, t) = [C{x, z) + A(z)] P{x, z, t) , (2) 

where C{x,z) = ^[U'{x) -\-V'{x)z] + q-§^ is the 
Fokker-Planck (FP) operator. The free evolution of 
the barrier noise is governed by the operator A(z) ~ 
1/t [{d/dz)z + Qd^/dz"^] for OUN or by the matrix A = 
(—7, 7; 7, —7) for DN. Initially the particle is located at 
the bottom Xb of the well and the quantity of interest 
is the mean first passage time (MFPT) through a given 
threshold Xthr located either at the top Xt or far from it 
on the other side of the barrier. 

A typical scenario of an escape event consists of two 
stages. For a long time t\j the particle fluctuates in the 
vicinity of the bottom of the well, being subjected to 
small random impacts of ^(t). If a large enough outburst 
of ^(t) occurs the particle will eventually surmount the 
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barrier almost immediately, during a short time U- The 
time- variation of the potential exerts only a negligible ef- 
fect on the first stage, but it can essentially modify the 
dynamics during the second one, when the particle inter- 
acts with the whole slope of the barrier. Any realization 
of ^{t), which has been supposed to bring the particle 
over the top of a static barrier, may turn out to be insuf- 
ficient if the barrier rises during the climbing stage. On 
the contrary, if the barrier decreases the particle does 
cross to the other side, but some smaller outbursts of 
£,{t) would also result in a successful escape. Because the 
rate of variation of the barrier shape depends on the cor- 
relation time of z{t), the relationshi p b etween tt and r 
appears to be crucial in the analysis |17l |. 

This discussion leads us to the central idea of the 
present approach - splitting the barrier noise into two 
independent components: 



Z{t) ^ Zg{t) + Zf{t) . 



(3) 



The slow one Zs is defined as the mean value of z over the 
time interval of climbing {t,t + tt) and over its possible 
realizations (marked by (...)) 



Zs{t) - ( ^ 



t+tt 



1 



dszis))^-{l-e-^)z{t), (4) 



where A = tt/r. It is supposed to be constant during the 
climbing stage, while its random character arises from the 
randomness of z(t). Hence Zg is governed by the same 
statistics as z but with the variance 



Dg 



tt A 



(1 



(5) 



Next, assuming that fast part Zf{t), which gives rapid 
fluctuations around Zg{t), can be treated as uncorrelated, 
one calculates its intensity Q / : 



Qf = Q 



1 
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(6) 



If z{t) is Gaussian it can always be written as the sum of 
two independent Gaussian components So, in OUN 
case both Zs{t) and Zf{t) are OUN's with correlation time 
r and they differ only in the form of their intensities 
(variances) Di — Qi/r {i — f,s). If r ^ one has 
Qf = Q, while the leading-order term of Qg reads Q/A^ 
so, for any a, it vanishes at least linearly with r. Thus 
one is left with only the fast part of z{t). In the opposite 
limit r ^ oo the leading term oi Qf becomes QA^/3, so 
Df vanishes at least linearly with 1/r, while Dg — D. 
Only the slow part of z{t) survives. One can check that, 
ignoring the dependence of Q on r, the intensities Q f and 
Qs are monotonic functions of r. While for r = one has 
the white-noise limit of z{t) with rapid fluctuations Zf{t), 
an increase of r increases the role of Zg at the expense of 
decrease of the intensity of zy, eliminating it completely 



as r — > cx). Thus the fundamental difference between 
Zg{t) and Zf{t) consists in the different regime of values 
of T in which they exist: Zg{t) occurs for t > tt, and 
hence fluctuates slowly, while Zf{t) persists for t < tt 
and varies rapidly. Only for t ~ it do they coexist. 

A similar summation property to that for Gaussian 
noise does not apply to the dichotomic noise — one can- 
not display a given dichotomic noise as the sum of two 
independent dichotomic noises. However, the great simi- 
larity between the statistical properties of OUN and DN 
suggests treating the DN case in the same way. The def- 
initions and Q involve the asymmetric character of 
two-state noise Zf(t) and its dependence on Zg{t), but 
for simplicity we assume, that both Zg(t) and Zf{t) are 
symmetric, independent dichotomic noises of zero mean. 
Since OUN and DN have the same correlation function 
the formulas (|4I6|I apply to the DN well. 

We should also determine the value of integration inter- 
val tt- For an unperturbed potential it equals the relax- 
ation time tr from the top to the bottom of the well, but 
fluctuations of the potential lead to far-from-equilibrium 
conditions, so that this equality does not hold |l9|. How- 
ever, we do not intend here to consider the relationship 
between the processes of climbing up and relaxing down 
the fluctuating barrier. Rather, we need a tool for cal- 
culating the order of the duration of the second stage of 
the escape event. It is enough to take for it the value 
of tr for a static barrier, which may be calculated as the 
MFPT from the top xt to the bottom xt of the well. It 
is shown in Fig. ^that our results depend almost unno- 
ticeably on the variation of tt within the range of tens 
of percent. A more careful analysis would require us to 
take into account, not only the mean value, but also the 
statistical distribution of relaxation times |15(|. 

Using the decomposition © the escape problem may 
be considered as a three-dimensional markovian process. 
Its joint probability distribution P{x, Zf,Zg,t) evolves ac- 
cordingly to the FP equation similar to Q but with 
two A's operators for Zf and Zg (with Qf or Qg instead 
of Q, respectively), and z = Zf + Zg in C{x,z). Such 
a formulation allows for a clear separation of different 
time-scales of the system dynamics. Since, by defini- 
tion, Zg remains constant while the particle climbs the 
barrier, its dynamics may be analyzed by the kinetic 
approach. On the contrary, Zf vanishes for r slightly 
greater than tr, but still for t ^ T, so rate theory ap- 
plies. Thus we seek the probability distribution in the 
form P{x, Zf, Zg,t) = p{x, Zf,t; Zg)p(zg,t) 0|. 

The fast equilibration process is described by the evo- 
lution oip{x, Zf,t;zg), which is governed by the equation 



d 

—p{x,zf,t;zg) = [C{x,Zf;Zg) + K{zf)]p{x, Zf ,t; Zg) , 

(7) 

where C{x, Zf;Zg) = ^ \U'{x; Zg) + V'{x)zf] + q-§^ and 
the slow component of barrier fiuctuations gives rise to 
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different forms of potential configurations U{x]Zs) = 
U{x) + V{x)zs. Following a standard method one looks 
for the quasi-potential $ being the dominant exponential 
term of the reduced (quasi-)stationary probability distri- 
bution {p{x,Zf;Zs))^^ = p{x\Zs) ^ exp Zs)/^] of 
For the DN case we obtain an equation 

[U'{x;zs)-'^'{x-zs)f 

^^[G{x)^\x-z,)~U\x-z,)] , (8) 
r 

whose middle (of the three always real) solution gives 
the quasi-potential. This equation is formally similar to 
the result of Reimann and Elston i who consider the 
case T ^T, however. The only difference is the form of 
diffusion function G{x) = 1 + {Qf/q)V'{x)^. In [l^ the 
total noise intensity Q is used, what gives an improper 
limiting value of $' for r — > cx) for a = 1. Here G{x) 
depends on Q^, which vanishes for any a as r — > cxd, so 
one obtains the exact expression $'(x; Zg) — > U'{x; Zg). In 
the opposite limit of r 0, the solution of © converges 
to the exact form U'{x)/G{x). This suggests to deal not 
with the quasi-potential but rather with an effective one 




FIG. 1: Relative mean escape time T/Ts versus r for DN case 
with a triangle barrier U{x) = 10(1 — |x|) and V{x) = 1 — \x\ 
confined to the interval ( — 1, 1) for Qo = 1, q = 1, xthr ~ 0, 
tr = 0.09 and a = 0, 0.25, 0.50, 0.75, and 1.0 from the bottom 
to the top on the left-hand-side, respectively. Solid lines show 
the exact results and the others the approximation (1131 for 
few values of tt- 



U'^ff{x] Zs 



<^'{x-Zs)G{x) 



(9) 



Finally, exploiting the well known form of the exact FP 
equation in the white-noise- limit |22|. one can write the 
effective FP operator 



d 



d 



d 



Ceffix;zs) = —U'^ff{x;z,) + q—^G{x) — ^G{x), 



dx 



dx 



(10) 

which governs the fast part of the evolution. 

A convenient way of finding the quasi-potential in the 
OUN case formulates the problem by means of path- 
integral or Hamiltonian techniques |23j. In general, the 
problem cannot be elaborated analytically, but asymp- 
totic expressions for small and large r are available 
[nl . To attempt an interpolation between the two 
limits of T we construct a 2-2 Pade approximant |i24i | 



-.- / / \ ix'. z.^) 

1 + r 2QfV'{xfW{x- Zs)/G{x) + W{x; ZsY 
l + T^Wix;zs)^/G{x) 



(11) 



with W{x;z,) = [G{x)/V'{x)][U'{x;Zs)V'{x)/G{x)]'. 
One can notice, that as a function of r the expres- 
sion ()ll|l has no singularities and monotonically increases 
wit h T. what is an anticipated property of quasi-potential 
hj,b3. As for DN, we may also introduce an effective 
potential ®. Using ^ calculation of the MFFT T{zs) 
for both types of barrier noise is straightforward. 

In the slow time scale the evolution of the system is 
governed by the Smoluchowski equation with a sink term 



—p(zs,t) = [A{zs) - k{zs)\ p{zs,t) 



(12) 



It describes stochastic switchings between the poten- 
tial configurations of different z^ and an escape process 
from each of them [k{zs) = ii/T{zs) with /i = 1/2 for 
Xthr — Xt, or II = 1 for xthr far from it]. One gets the 
mean escape time integrating p{zs,t) over t G (0, oo), 
and summing/integrating over Zg for DN/OUN. For the 
dichotomic switching the result is immediate: 



T = 



2T+T- + fiT{T+ +T-^) 



2fJ.T 



(13) 



where T± are the MFPT's for U±{x) = U{x) ± ^/lTsV[x), 
respectively. Although (|13|l resembles the well-known so- 
lution |3, 123 of a very simple set of equations l(T^ which 
constitutes the long-r approximation of the problem j 
the dependence of 7!|_ and T_ on Qfir) involves also the 
fast part of the dynamics in the formula H13|) . 

The problem is much more complicated in the OUN 
case. To the best of the author's knowledge there is no 
universal approximation of (|12|l valid for any r |2^ . One 
may calculate asymptotic expressions for small and large 
T 0, l^^l and construct a Pade approximant to inter- 
polate in between; however, the complicated exponential 
dependence of expansion terms on the amplitude of fiuc- 
tuations yields a very bad approximation 28]. Hence, in 
what follows, we solve (|12() numerically. 

To test the method we take the triangular barrier 
model |i] with DN. In Fig. [T] we plot T{t)/Ts {Ts is 
the MFPT for a static barrier) for the exact analytical 
results and for the present method, in each case for few 
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10"^ 10"' 10"' 10° 10^ 10' 10\ 10^ 

FIG. 2: Relative mean escape time versus r for OUN 

case and the system witli U{x) = a;'*/4 — V{x) — 

U{x) + 1/4 for \x\ < 1 and V{x) = elsewhere, for q = 0.08, 
Qo = 0.8, tt — tr/2 — 1.26, and three values of a. The lines 
present our approximation and markers are from the numer- 
ical simulation of Q. 

values of a. The relaxation time calculated from the ex- 
act formula for the MFPT from a;t = to Xf, = 1 
equals tr = 0.09. We show three sets of curves with 
tt = tr, tt = tr/2 and tt = tr/4:, respectively. The agree- 
ment with the exact plot is very good, but in the interval 
10"'^ < r < 10~^ our method gives slightly lower values. 
We have found the smallest deviation for tt = tr/2, but 
even when tt is twice larger or smaller the difference is 
still not very significant. This validates the way we esti- 
mate the interval of integration tt in Q). For simplicity 
in the next example we use tt = tr/2, but to be more pre- 
cise, for each system a careful analysis of its best value 
should be done Hi]. In Fig. El we display T{t)/Ts for 
OUN case and three values of a. The agreement between 
the theory and numerical simulation of is very good, 
but also with some underestimation in the region of the 
resonant activation minima. The results for other sys- 
tems and other values of parameters are also excellent 

m 

To conclude, we have presented a method of an inves- 
tigation of thermal activation in the presence of barrier 
fluctuations for arbitrary duration of their correlation. 
Dividing the barrier noise into two components - the slow 
and fast ones - we can separate two time scales of the 
evolution of the system for any value of r and use both 
rate and kinetic approaches in the analysis without any 
sewing procedure. The noise division is done through an 
averaging over a finite interval of time tt Q , hence we call 
the approach a partial noise-averaging method (PNAM). 
For a dichotomic perturbation the formula (|13|l together 
with the MFPT obtained for the FP operator (HHJ pro- 
vides for the first time the analytical expression for the 
dependence T{t) for any t S [0,cx)), for arbitrary po- 
tentials U{x) and V{x), and a large class of noises. For 



the OUN we have been obliged to use a computer at the 
final step, but the accordance of the present result with 
the full-numerical ones confirms the power of PNAM. Al- 
though the method is presented in terms of MFPT, it can 
be expressed by means of any of the standard approaches 
Q to the activation process. We hope also, that the pre- 
sented idea of splitting the noise could be useful in other 
problems where different time-scales coexist, making the 
proposed approach valuable for many applications. 
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